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摘 要 : 以 叶 尔 羌 河 流域 上 游 河源 区 为 研究 区 ,构建 说 入 了 冰川 模块 的 SWAT 分 布 式 水 文 模型 ,对 
卡 群 水 文 站 1968 一 2017 年 逐 月 径流 进行 模拟 ,评价 该 模型 在 研究 区 内 的 适用 性 ,基于 模拟 结果 分 


析 研 究 区 冰川 径流 年 际 变化 与 年 内 分 布 情况 ,并 定量 核算 气候 和 下 热 本 


变化 对 径流 变化 的 贡献 


率 。 结 果 表 明 :该 模型 在 研究 区 径流 模拟 中 具有 良好 的 适用 性 ,校准 期 (1968 一 1992 年 ) 与 验证 期 
(1993 一 2017 年 ) 的 决定 系数 ( 尼 ) 分 别 为 0.77 和 0.86, 纳 什 系数 (NSE ) 分 别 为 0.76 和 0.85, 均 方 根 误 
差 与 实测 值 标准 差 的 比值 (RSR ) 分 别 为 0.49 和 0.38, 偏 差 百 分 比 (PBIAS) 分 别 为 -7.4% 和 0.6%。 经 
模拟 ,研究 区 近 50 a 冰川 径流 量 总 体 呈 增加 趋势 , 且 其 占 总 径流 量 的 比重 约 为 51.1%; 年 内 冰川 产 


流 主要 发 生 在 6 一 9 月 , 占 全 年 冰川 径流 量 的 90.0% 以 上 。 气 候 和 下 垫 理 


变化 分 别 使 月 径流 量 增加 


6.62 m"-s 和 0.41 ms ,其 中 气候 变化 的 贡献 率 为 94.2%, 即 气候 变化 对 研究 区 径流 的 影响 占 主 时 
地 位 。 研 究 成 果 可 为 理解 研究 区 的 历史 径流 变化 成 因 及 预测 未 来 径流 的 演变 趋势 提供 科学 依据 。 
X 键 词 : 气候 变化 ， 下 垫 面 变化 ，SWAT 模 型 ;冰川 径流 模拟 ; 叶 尔 羌 河流 域 


文章 编号 : 


气候 和 下 垫 面 变化 是 影响 流域 水 文 循环 过 程 
改变 的 主要 因素 ,其 中 径流 演变 对 不 同时 空 尺 度 下 
气候 和 下 垫 面 变 化 的 响应 研究 是 水 资源 规划 及 其 
管理 的 基础 "。 目 前 在 该 研究 领域 主要 采用 的 研究 
方法 以 对 比分 析 法 ”统计 分 析 法 和 水 文 模型 法 “ 
为 主 。 其 中 ,基于 物理 过 程 模 拟 的 分 布 式 水 文 模型 


究 工具 之 一 *。 

近年 来 随 着 SWAT 模 型 的 发 展 日 趋 成 熟 , 针 对 
该 模型 的 应 用 也 在 不 断 完善 ,例如 ,Rahman 等 "在 
瑞士 罗 纳 河上 游 山区 应 用 SWAT 模 型 研究 了 径流 对 
区 域 气候 模式 输出 的 响应 ; 郭 军 庭 等 "在 潮 河流 域 
应 用 SWAT 模 型 进行 月 尺度 径流 模拟 ,量化 区 分 了 


能 够 充分 反映 流域 的 空间 异 质 性 ,较为 准确 地 刻画 
流域 水 文 过 程 , 并 揭示 水 文 循环 过 程 的 内 在 作用 机 
理 。 该 方法 已 广泛 应 用 于 气候 和 下 垫 面 变 化 及 两 
者 协同 变化 对 径流 影响 的 研究 中 “”"。SWAT 模 型 
(Soil and water assessment tool ) 作 为 一 种 基于 GIS 的 
半分 布 式 流域 水 文 模型 ,已 应 用 在 我 国 多 个 流域 的 
径流 模拟 研究 中 ,特别 是 在 实测 数据 较为 缺乏 的 西 
北 内 陆 河 流域 ,该 模型 也 是 径流 过 程 模拟 的 重要 研 
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下 热 面 和 气候 变化 对 径流 变化 的 影响 ; 赵 杰 等 "在 
乌鲁木齐 河上 游 汇 水 区 应 用 SWAT 模 型 进行 日 径流 
模拟 ,发 现 模拟 结果 与 观测 流量 过 程 线 拟 合 程度 较 
好 ; 刘 全 刘 等 “和 祁 敏 等 在 水 河 和 阿克苏 河流 域 
分 别 结合 GIS 和 加 入 融雪 模块 的 SWAT 模 型 对 多 年 
月 径流 进行 了 模拟 ,发 现 结合 GIS 和 融雪 模块 后 的 
模型 模拟 精确 度 更 高 。 但 目前 国内 外 学 者 多 是 从 
降雨 或 融雪 的 角度 对 径流 进行 的 模拟 研究 ,对 于 冰 
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川 水 文 过 程 模拟 的 相关 研究 涉及 较 少 ,例如 在 塔 里 
木 河流 域 源流 区 , 张 一 驰 等 的 和 黄 粤 于 0 等 分 别 应 
用 MIKE SHE 模型 对 开 都 河流 域 径流 进行 了 模拟 研 
究 ; 王 月 等 ”在 阿克苏 河流 域 应 用 含 融雪 结构 的 新 
安江 模型 对 出 山口 径流 进行 模拟 ,模拟 结果 精度 有 
待 提 高 ; 罗 敏 等 应 用 MIKE SHE 模型 对 未 来 气候 
情景 下 的 和 田 河 流域 日 径流 过 程 进行 了 模拟 , 率 定 
及 验证 期 的 模拟 结果 均 在 可 接受 范围 内 。 此 外 ,在 
干旱 内 陆 河 流域 ,许多 以 冰雪 融 水 为 主要 补给 源 的 
河川 径流 模拟 结果 也 不 够 理想 ,这 些 人 研究 成 果 无 法 
为 流域 管理 者 合理 制定 水 资源 管理 措施 提供 决策 
XE, 

本 文 以 叶 尔 芜 河流 域 上 游 河 源 区 为 研究 区 , 利 
用 冰川 编目 数据 修正 研究 区 内 的 冰川 空间 分 布 情 
况 ,构建 能 够 模拟 冰川 水 文 过 程 的 SWAT 模 型 对 月 
径流 过 程 进行 模拟 ,分 析 气 候 和 下 垫 面 变化 对 径流 
演变 的 贡献 率 ,以 期 为 叶 尔 芜 河流 域 上 游 河 源 区 的 
生态 保护 水 资源 管理 及 其 相关 政策 措施 的 制定 提 
供 科学 依据 。 


1 研究 区 数据 与 方法 


1.1 研究 区 概况 

叶 和 尔 芜 河流 域 上 游 河 源 区 地 处 喀 喇 昆仑 山北 
Jit Ho BA bs 74927! —78924' E, 35?26' -38?18' N) , 流 
域 面 积 为 46404.87 km’, 占 全 流域 总 面积 的 46.9%， 
人 研究 区 地 势 南 高 北 低 , 以 高 山 积 雪 带 、 亚 高 山 、 高 山 
草 旬 带 和 低 山 丘陵 带 等 地 带 为 主 20 ;海拔 在 1447~ 
8537 之 间 , 具 有 明显 的 大 陆 气 候 特 征 ( 图 1)。 研 
究 区 境内 共 发 育 有 第 四 纪 冰 川 3059 条 ,冰川 总 面积 
5315 km , 总 储量 达 684.5 km ,年 冰川 融 水 量 达 
38.24x10* ms201。 研 究 区 多 年 平均 气温 为 -10.2 C, 
当地 降水 主要 受 地 形 影 响 ,高 山 带 年 降水 量 约 500~ 
700 mm, 中 高 山 带 年 降水 量 约 300 mm, 中 低 山 带 年 
降水 量 约 130 mm。 人 研究 区 水 系 属 典型 冰川 融雪 补 
给 型 河流 ,其 中 干流 长 约 585 km ,年 径流 量 66x108 
m ,其 中 冰雪 融 水 径流 占 比 77.4%2。 据 研究 区 下 
游 出 山口 处 水 文 站 一 一 卡 群 站 (地 理 坐 标 76"54' 工 ， 
37°59' N; IR 1450 my) 的 实测 资料 分 析 显 示 ,该 地 
区 径流 量 年 内 分 布 极 不 均匀 ,其 中 6 一 8 月 径流 量 约 
占 年 径流 量 的 73.8%。 径 流出 山 后 ,经 叶 尔 羌 河 流 
域 中 下 游 并 最 终 与 和 和 田 河 一 起 汇 人 塔里木 河 干流 。 


>z 


^ 卡 群 水 文 站 点 
。 0.5? 虚拟 气象 站 点 
a 塔什库尔干 气象 站 点 


图 1 叶 尔 羌 河流 域 上 游 河源 区 概况 
Fig. 1 Overview ofthe headwaters of the upper Yarkant 


River Basin 


12 基础 数据 
1.2.1 DEM 数据 及 冰川 编 卓 数据 数字 高 程 模型 
(DEM) 数 据 来 源 于 国际 科学 数据 服务 平台 ,分 辨 率 
为 90 m, 主 要 用 于 提取 地 形 信息 .计算 水 流 方向 和 
集 水 面积 .提取 河 网 、 划 分 子 流域 。 

研究 区 第 一 次 冰川 编目 数据 来 源 于 1962 一 
1977 年 航 测 地 形 图 ,第 二 次 冰川 编目 数据 来 源 于 
2004—2011 年 Landsat TMETM+ 和 ASTER 遥感 影 
像 。2 次 冰川 编目 数据 分 别 用 于 确定 校准 期 和 验证 
期 内 的 冰川 空间 分 布 情况 。 
1.4.2 下 垄 面 数 据 人 研究 区 下 垫 面 数据 包括 我 国境 
内 和 境外 2 个 部 分 。 国 内 部 分 以 地 理 空间 数据 云 
(http://www.gscloud.cn/)1990、2015 年 Landsat 数 据 作 
为 基础 数据 源 ,分 辩 率 为 30 m。 境 外 部 分 数据 以 欧 
盟 联合 研究 中 心 的 全 球 土地 覆盖 数据 集 (CLC2000) 
作为 数据 源 Chttps://forobs.jrc.ec.europa.eu/products/ 
glc2000/products.php) 。 本 文 利用 第 一 次 和 第 二 次 
的 冰川 矢量 数据 分 别 蔡 换 1990 年 和 2015 年 下 垫 面 
类 型 中 的 冰川 类 型 数据 (CNJX)。 依 据 SWAT 模 型 
标准 建立 下 垫 面 类 型 索引 表 并 进行 重 分 类 ( 表 1)， 
最 终 得 到 模型 所 需 的 下 热 面 类 型 图 (图 2)。 
1.2.3 土壤 数据 研究 区 国内 部 分 的 土壤 数据 以 中 
国 科 学 院 南 京 土壤 所 发 布 的 1:1000000 的 土壤 数据 
库 为 数据 源 ,境外 部 分 数据 以 联合 国 粮农 组 织 
(FAO) 和 维也纳 国际 应 用 系统 研究 所 (IIASA) 构 建 
的 世界 土壤 数据 集 (Harmonized world soil database, 
HWSD ) 为 数据 源 (http://westdc.westgis.ac.cn)。 本 文 
利用 第 一 二 次 冰川 矢量 数据 替换 了 原 数 据 中 的 冰 
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表 1 下 扑面 重 分 类 名 称 及 面积 占 比 


Tab.1 Names and area proportions of underlying surface 


川 雪 被 类 型 (ICE)。 各 项 土壤 属性 参数 通过 SPAW 
软件 及 检索 中 国土 壤 数 据 库 (http://vdb3.soil.csdb. 


ee ae cn) 得 到 ,依据 SWAT 模 型 标准 对 原 数 据 进行 重 分 类 
I pm REM (62) ,最 终 得 到 模型 所 需 的 土壤 类 型 图 (图 3)。 
Med i 124 气象 水 文 数据 ”由 于 研究 区 内 仅 有 一 个 国家 
P i ien "n 气象 站 (塔什库尔干 站 ), 仅 用 该 站 点 的 实测 数据 无 
31 ” ”高柳 羡 草 地 XJCD 9.8 法 反映 整个 区 域 的 气象 环境 。 因 此 ,本 文 的 气温 和 
32. quU XJCY 17.5 降水 数据 采用 中 国 地 面 降水 和 气温 日 值 0.5° 格 点 数 
33 mmus WHM is 据 集 (V2.0) 资 料 ,来 源 为 中 国 气象 数据 共享 网 (http: 
41 dp B XJHQ 01 //data.cma.cn)。 该 数据 集 是 由 我 国境 内 751 个 观测 
Wi o CE b 让 的 观测 数据 插值 形成 =” al ETE SUI ERIS A 
PU $5 析 检 查 ,可 以 较 好 地 反映 研究 区 的 气温 和 降水 空间 
62 XE XJGB 16 分 布 特征 。 本 文 用 于 模型 率 定 与 验证 的 径流 数据 
66 WATE XJLY 48.7 来 源 于 卡 群 水 文 站 1968 一 2017 年 实测 月 径流 资料 。 
(a) 1990 年 (b) 2015 年 
下 垫 面 类 型 
D e PY ll 中 禾 盖 度 草 地 。 pum mm 冰川 E 戈壁 
mmn EH 高 覆盖 度 草地 加 | 低 覆 盖 度 草地 — 水 库 坑 塘 — 国 四 城镇 用 地 。 Emp 裸 岩 地 
图 2 叶 尔 羌 河 流域 上 游 河 源 区 下 垫 面 分 布 
Fig.2 Distribution of underlying surface types in the headwaters of the upper Yarkant River Basin 
表 2 土壤 重 分 类 名 称 及 面积 占 比 
Tab.2 Names and area proportions of soil reclassification 
亚 类 名 称 土 类 名 称 土壤 亚 类 代码 重 分 类 名 称 面积 占 比 /% 
RRE RE 23120171 HDTXJ 41.2 
标 钙 土 棕 钙 土 23113101 ZGTXJ 9.2 
淡 标 钙 土 标 钙 土 23113102 DGTXJ 54 
HRE RE 23114121 ZMTXJ 3.3 
含 盐 石 质 土 石 质 土 23115185 HYSZTXJ 1.4 
Tg Tg 23120102 CZTXJ 0.9 
寒 钙 土 寒 钙 土 23120122 HGTXJ 2.0 
冷 钙 土 冷 钙 土 23120132 LGTXJ 18.7 
寒 漠 土 RRE 23120151 HMTXJ 5.3 
冰川 雪 被 冰川 雪 被 23127101 ICE 12.6 
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(a) 1990 年 


0 50km 
LL 


(b) 2015 年 


0 50km 
| 


土壤 类 型 
ES 棕 钙 土 mm dx Eu t At mm GUI 
mm 淡 棕 钙 土 Em 含 盐 石 质 士 BSE EAE 冰川 雪 被 


图 3 叶 尔 羌 河 流域 上 游 河源 区 土壤 类 型 分 布 
Fig. 3 Distribution of soil types in the headwaters of the upper Yarkant River Basin 


1.3 研究 方法 
1.3.1 降水 、 气 温和 径流 变化 趋 芬 分 析 采用 一 元 
线性 回归 法 后 对 研究 区 气候 特征 及 年 径流 历史 变 
化 趋势 进行 分 析 ,并 采用 非 参数 Mann-Kendall 法 对 
其 变化 趋势 进行 显著 性 检验 。 
1.3.2 冰川 模块 算法 

Mz G ERRAT anic) > T> Tar 


m 


0 ? T< T. aiciss q 


式 中 ;MW 为 冰川 每 日 消融 量 (m) ; Fu Zr vi i f D] 
子 ;Ri 为 冰 辐 射 融 化 因子 ;1, 为 潜在 太阳 直接 辐射 ; 
+Rilw 为 考虑 辐射 因素 影响 的 冰川 消融 因子 ;7 为 
日 均 温 (% ) ; T, LL 为 发 生 冰 川 消融 的 阔 值 温度 
(C); 
V2cS" (2) 

式 中 :V 为 冰川 体积 (mw) ic 为 常数 ;$ 为 冰川 表面 积 
(m);7 为 无 量 纲 标 度 系 数 。 

F=W. xp + sin 2 e : sp] (3) 
式 中 :有 为 冰川 物质 累积 量 ; 取 为 积 雪 水 当量 ; B 为 
基础 累积 系数 ;为 某 天 的 序号 。 
1.3.3 水 文 过 程 模拟 本文 基 于 SWAT 模 型 对 叶 尔 
羌 河 流域 上 游 河源 区 的 径流 进行 模拟 ,进而 量化 分 
析 气 候 和 下 垫 面 变化 对 当地 径流 的 影响 。SWAT 模 
型 是 由 美国 农业 部 农业 研究 中 心 (USDA-ARS ) 开 发 
的 1 个 基于 物理 过 程 的 半分 布 式 水 文 模型 ,目前 已 
被 广泛 应 用 于 流域 水 文 模拟 流域 管理 及 非 点 源 污 
染 控 制 等 方面 ”1。SWAT 模 型 建 模 过 程 大 致 包含 
3 个 部 分 :首先 基于 DEM 数 据 对 研究 流域 进行 流域 
离散 化 ,然后 基于 模型 中 输入 的 土地 利用 、 土 壤 数 


据 及 坡度 等 进行 水 文 响应 单元 (HRUs) 划 分 ,最 后 根 
据 输入 的 气象 数据 模拟 流域 的 水 文 过 程 ,并 根据 模 
拟 结果 对 模型 参数 进行 调整 和 验证 ,人 研究 流程 如 下 
图 4 所 示 。 本 文 将 Ying 等 ”开发 的 冰川 径流 的 冰川 
模块 馈 入 SWAT 模 型 ,对 研究 区 冰川 融 水 径流 进行 
模拟 。 

研究 区 以 200 km 为 集 水 冰 值 ,将 流域 划分 为 
112 个 子 流域 ,2379 个 水 文 响应 单元 (HRUs) ,模拟 
近 50 a(1968 一 2017 年 ) 逐 月 径流 ,其 中 1968 一 1992 
年 为 模型 校准 期 ,1993 一 2017 年 为 模型 验证 期 。 采 
用 SUFI-2 不 确定 性 分 析 算 法 (Sequential uncertainty 
fitting-version 2) 对 模型 参数 进行 敏感 性 分 析 ,采用 
自动 率 定 和 手动 校准 相 结合 的 方式 对 参数 进行 率 
定 。 模 拟 结果 的 评价 标准 选用 决定 系数 CR) .纳什 
系数 (Nash-Sutcliffe efficiency, NSE) 、 均 方 根 误差 与 
实测 值 标准 差 的 比值 (RMSE-observations standard 
deviation ratio, RSR) 和 偏差 百分比 (Percent bias, 
PBIAS)4 项 指标 ,定义 公式 如 下 : 


Po E = ez) 


i- 


(4) 


R 


(5) 
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任 才 等 : 气候 与 下 垫 面 变 化 对 叶 尔 羌 河 源流 径流 的 影响 


I3 f DEM 数 据 | Aat qM 
i 水 ; | | Ju | 区 分 ! 
i Oow) | 并 量 | 
Aj i ; NECS : i 
PINK i ENE ILALL 
T i de 一 > 型 一 一 一 > RD i 
| x EM 土壤 数据 Ha a "s 
1 : i 1 1 i yn| ' 
i 化 | EN | 源 径 ， 
rt 气象 水 文 数据 I E 
ES (d i ; 影响 
i 分， BAS 0 0 d — 0 odo 0 Xu 
er i 冰川 模块 j 
! SWAT 模 型 建立 及 本 土 化 -一 
图 4 研究 流程 图 
Fig. 4 Research process 
*( Om Q;") 变化 对 径流 的 影响 程度 ,计算 公式 如 下 : 
PBIAS = 一- — — — x100 (7) AO 2 AQC* AQ, (8) 
obs A 
20 B. Oc y 1999; (9) 
i "^ | s He Hen AAA SNL E) C+L 
式 中 : 097" 为 时 间 步 长 i 内 模型 模拟 的 河道 径流 量 AQ 
— L 
(m) x QP 为 时 间 步 长 ;内 水 文 站 实际 观测 的 河道 径 


WEE (m?) ; 0 IL Qo 分 别 为 模拟 时 间 段 内 模拟 和 


式 中 : AO 为 径流 在 气候 和 下 垫 面 变化 共同 作用 


观测 的 平均 径流 量 (m); STDEV,,. 为 实测 值 标准 
差 , 反 映 实 测 值 相 对 于 平均 值 的 离散 程度 。 

NSE 系数 用 来 评价 模型 模拟 的 精度 , 取 值 范围 
为 -w~1, 越 接近 1, 说 明 模 型 模拟 效果 越 好 ,一 般 认 
为 NSE>0.75 ,模拟 结果 为 “ 优 ”。RSR 指标 越 接 近 于 
0 表示 模拟 效果 越 好 ,一 般 认 为 RSR<0.5 ,模拟 结 
为 “ 优 ”。PBIAS 表示 模拟 的 平均 值 高 于 或 低 于 实 
测 值 , 越 接近 0 表明 模拟 效果 越 好 ,一 般 认为 
PBIAS<+10% ,模拟 结果 为 " 优 "。 尼 取 值 范围 为 0~ 
1, 表 明 模 拟 值 与 实测 值 之 间 的 线性 相关 程度 和 相关 的 
方向 , 越 接 近 1, 表 明 模 拟 值 与 实测 值 相关 性 越 好 ”。 
1.3.4 气候 和 下 垫 面 变 化 及 二 者 协同 作用 对 径流 的 
影响 ”将 研究 期 内 的 气候 和 下 扑面 数据 相互 交叉， 
并 设计 出 SS S 、S 共 4 种 子 情景 ( 表 3) ,分析 对 比 
不 同情 景 下 的 径流 模拟 结果 ,以 量化 气候 和 下 热 面 


FAS 25 E E T (Wm? s) ; AQ, 为 径流 在 气候 变化 作 
用 下 的 变化 量 (ms ; AQ, 为 径流 在 下 垫 面 变化 
作用 下 的 变化 量 (m…s'); Be 和 B, 为 气候 变化 和 下 
垫 面 变化 对 径流 影响 的 贡献 率 (% )。 


2 结果 与 分 析 


21 降水 量 、 气 温和 径流 变化 特征 

根据 塔什库尔干 气象 台 站 及 卡 群 水 文 站 相关 
实测 数据 ,对 研究 区 1968 一 2017 年 降水 量 .气温 和 
径流 序列 历史 变化 趋势 分 析 和 显著 性 检验 结果 (图 
5 图 6 和 图 7) 表 明 , 近 50a 研 究 区 年 均 降 水 量 增加 
了 39.8 mm, 年 增长 率 为 1.0% ;年 均 气 温 上 升 了 
2.44 ,年 增长 率 为 1.2% ;年 径流 量 增加 了 17.20x 
10 mi ,年 均 增 长 率 0.5% ;由 Mann-Kendall 趋势 检验 


表 3 气候 和 下 垫 面 变化 对 径流 影响 的 情景 设置 


Tab.3 Scenario settings of the influences of climate and underlying surface change on runoff 


情景 气候 数据 年 份 下 热 面 数据 年 份 模拟 径流 结果 
S, (基础 情景 ) 1966 一 1992 1990 Qo 
S, (气候 变化 情景 ) 1993 一 2017 1990 Q. 
S, (下 热 面 变化 情景 ) 1966 一 1992 2015 Q. 
S, (协同 作用 情景 ) 1993 一 2017 2015 Q: 
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年 降水 量 /mm 


05 LU 1. 1. 1 1. 1 1 1 ] 
1968 1973 1978 1983 1988 1993 1998 2003 2008 2013 2018 
图 5 塔什库尔干 站 年 降水 量变 化 趋势 


Fig. 5 Change trend of annual precipitation in 


Taxkorgan Station 


法 知 ,在 w=0.05 的 显著 性 水 平 下 ,降水 量 .气温 和 径 


流 检验 均 通 过 信和 度 99% 显 著 性 检验 ,可 见 三 者 的 增 E 
加 趋势 均 显 著 。 E 
22 模型 参数 选取 与 模拟 结果 分 析 ko 
22. 模型 参数 选取 与 率 定 ”运用 SUFI-2 算 法 对 相 ü 


关 模 型 参数 进行 敏感 性 分 析 , 从 中 选取 对 于 径流 模 
拟 精度 影响 较 大 的 13 个 参数 进行 自动 率 定 ,结合 研 
究 区 的 冰川 水 文 过 程 并 参考 Zhang 等 ”的 人 研究 成 
T ,采用 手动 校准 法 再 次 校准 参数 ,直到 模拟 结 
均 符 合 各 项 评价 指标 。 率 定 的 具体 参数 及 最 优 值 
见 表 4。 

2.2.2 模拟 结果 评价 与 分 析 ”本文 选取 1966 一 1967 
年 为 模型 模拟 的 预 热 期 ,未 输出 模拟 结果 。 以 
1968 一 1992 年 为 模型 校准 期 ,1993 一 2017 年 为 模型 
验证 期 ,对 研究 区 的 月 径流 进行 模拟 。 模 型 在 校准 


一 - 日 最 高 气温 -- 日 最 低 气温 


1968 1973 1978 1983 1988 1993 1998 2003 2008 2013 2018 


图 6 塔什库尔干 站 日 最 高 .最低 气温 变化 趋势 


Fig. 6 Variation trend of daily maximum and minimum 


temperature in Taxkorgan Station 


40 
1968 1973 1978 1983 1988 1993 1998 2003 2008 2013 2018 
年 份 
图 7 卡 群 站 年 径流 变化 趋势 
Fig.7 Change trend of annual runoff variation in Kaqun 


Station 


期 和 验证 期 的 月 值 模拟 结果 如 图 8 和 图 9, 参 数 化 评 
价 指 标 结果 如 表 5 所 示 。 

图 8 和 图 9 表明 ,模型 在 加 入 冰川 模块 后 ,校准 
期 和 验证 期 径流 模拟 曲线 和 实测 曲线 能 较 好 地 拟 


表 4 率 定 参数 及 其 最 优 值 


Tab.4 Calibration parameters and their optimal values 


参数 模块 参数 名 称 参数 意义 参数 范围 参数 最 优 值 

径流 ESCO 土壤 蒸发 补偿 系数 0.01, 1.0] 0.6 
ALPHA_BF 基 流 系数 0.1, 1.0] 0.036 
CH_K2 主 河道 河床 有 效 水 力 传导 度 0, 50] 6.5 
GWQMIN "diio" zz B 1, 200] 100 
GW. DELAY 也 下 水 延迟 时 间 0, 180] 60 

冰川 Fy 融 冰 系数 2.6, 13] 5.5 
Rie 冰 辐 射 系数 0, 0.1] 14.4x10? 
bBo 基础 蓄积 系数 0.001, 0.006] 0.003 

积 雪 SMFMX 6 月 21 日 融雪 度 日 因子 0, 10] 6.5 
SMFMN 12 月 21 日 融雪 度 日 因子 0, 10] 1.0 
SFTMP 降雪 阔 值 温度 -5. 45] 1.0 
SMTMP TRE E (Et JE -5, +5] 0.5 
TIMP JRE BER Jes RC 0.1, 2.0] 1.0 
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E FE: 气候 与 下 垫 面 变 化 对 叶 尔 羌 河 源流 径流 的 影响 
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图 8 卡 群 站 月 平均 流量 模拟 结果 (校准 期 ) 


Fig. 8 Simulation results of monthly average flow in Kaqun 


Station (calibration period) 
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图 9 卡 群 站 月 平均 流量 模拟 结果 (验证 期 ) 
Fig.9 Simulation results of monthly average flow in Kaqun 


Station (verification period) 


AEH , RRA UK THES ERES S SW AT PRALA 
效果 和 良好, 能够 较为 准确 地 刻画 研究 区 的 多 年 月 径 
流 过 程 。 人 研究 区 个 别 年 份 (1971 1973 .1977 1988 
年 和 1994 年 ) 的 冰 湖 溃 决 洪水 守 导 致 出 山口 汛期 的 
径流 模拟 结果 误差 相对 较 大 ,其 它 年 份 的 流量 峰值 
与 基 流 的 模拟 结果 均 较 为 准确 ,表明 模型 的 结构 较 
为 可 靠 ,模型 相应 参数 选取 合理 ,可 依据 模型 的 输 
出 结果 作 进 一 步 分 析 研究 。 

表 5 为 模型 校准 期 .验证 期 和 全 时 段 内 的 月 径 


流 统计 评估 结果 ,其 中 校准 期 月 径流 模拟 值 为 
218.10 m°- s”, EKMA 14.96 m^ s"! ,模拟 结果 相 
对 偏 大 ;验证 期 月 径流 模拟 值 为 224.72 m^ s" , 比 实 
测 值 低 1.34 mx…s ,表明 模型 在 验证 期 的 模拟 结 
与 实际 情况 更 相符 ;而 全 时 段 月 径流 模拟 值 为 
221.41 mms, 比 实测 值 略 高 6.81 ms ,但 全 时 段 的 
模拟 结果 能 基本 反映 实际 径流 的 变化 趋势 。 从 各 
项 评价 指标 来 看 ,校准 期 和 验证 期 的 评价 结果 均 为 
“ 优 ”, 其 中 校准 期 分 别 为 :R=0.77、NSE=0.76、RSR= 
0.49 .PBIAS=-7.37% ,表明 模拟 值 与 实测 值 的 拟 合 
程度 较 好 ;验证 期 分 别 为 :R=0.86、NSE=0.85、RSR= 
0.38 .PBIAS=0.59% ,表明 验证 期 模拟 结果 更 优 于 校 
准 期 。 
2.2.3 冰川 水 文 过 程 分 析 依据 分 布 式 水 文 模拟 模 
型 输出 结果 对 其 冰川 水 文 过 程 进 行 年 际 变化 趋势 
分 析 ( 图 10) ,结果 表明 : 近 50 a 研究 区 模拟 年 径流 
量 总 体 呈 缓慢 上 升 趋势 ,其 中 最 大 和 最 小 年 径流 量 
分 别 为 2012 年 的 99.55x10: mm 和 1996 年 的 48.59x 
10* m ,多 年 平均 年 径流 量 68.26x10: mm 。 而 冰川 径 
流 与 年 径流 模拟 过 程 的 变化 趋势 基本 一 致 , 呈 上 升 
趋势 。 年 最 大 和 最 小 冰川 径流 量 分 别 为 44.99x10 
m'(2001 4E.) fI 25.0x10* ms(1972 年 ) ,多 年 平均 冰川 
Í% fi E N 34.87 X 10* m, 占 多 年 平均 总 径流 量 的 
51.1%。 近 50a 冰 川 径流 量 总 体 呈 增长 的 趋势 ,年 均 
增长 率 约 为 0.2% ,其 占 总 径流 量 的 比重 在 34.0%~ 
70.4% 之 间 。 

研究 区 实测 径流 .模拟 径流 和 冰川 径流 年 内 分 
布 模拟 (图 11) 结 果 表 明 :研究 区 月 冰川 径流 模拟 结 
果 占 月 径流 模拟 总 量 的 比重 约 为 11.9%~74.0% ,其 
中 2 月 占 比 最 小 ,8 月 占 比 最 大 。 在 2 月 到 8 月 间 , 随 
着 气温 的 上 升 ,冰川 径流 量 呈 逐 月 增加 的 趋势 ;从 8 
月 至 翌年 2 月 ,冰川 径流 量 呈 逐 月 减少 的 趋势 。 在 
12 月 到 翌年 4 月 间 ,研究 区 的 冰川 径流 量 为 1.59x 
10° mm , 仅 占 年 冰川 径流 量 的 2.8%。 研 究 区 的 冰川 
产 流 主要 发 生 在 夏季 的 6 一 9 月 , 占 全 年 冰川 径流 的 


m5 卡 群 站 月 平均 径流 模拟 结果 评价 


Tab.5 Evaluation of monthly runoff simulation results in Kaqun Station 


时 段 实测 径流 均值 /as 模拟 径流 均值 ns” R NSE RSR PBIAS/% 
校准 期 (1968-01 一 1992-12) 203.14 218.10 0.77 0.76 0.49 -7.37 
验证 期 (1993-01 一 2017-12) 226.06 224.72 0.86 0.85 0.38 0.59 
全 时 段 (1968-01 一 2017-12) 214.60 221.41 0.82 0.81 0.44 -3.17 


注 : 尼 为 决定 性 系数 ,NSE 为 纳什 系数 ,RSR 为 均 方 根 误差 与 实测 值 标 准 差 的 比值 ,PBIAS 为 偏差 百分比 。 
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图 10 模拟 年 径流 ,模拟 冰川 径流 、 冰 川 径流 贡献 率 年 际 变 
化 趋势 
Fig. 10 Interannual variation trend of simulated annual 
runoff, simulated glacier runoff and contribution rate of 


glacier runoff 
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图 11 实测 径流 ,模拟 径流 ,冰川 径 流 和 冰川 径流 贡献 率 年 
内 分 布 特征 


Fig. 11 Distribution characteristics of the measured runoff, 
simulated runoff, glacier runoff and 


glacier runoff contribution during the year 


91.1%, 其 中 仅 8 月 的 比重 高 达 41.2%。 
2.3 气候 和 下 垫 面 变化 对 径流 演变 的 影响 

由 气候 及 下 垫 面 变 化 对 径流 影响 的 定量 分 析 
结果 可 知 ( 表 6) :基础 情景 (S50) 模拟 月 径流 Qo W 
217.67 m? s! ,气候 变化 情景 (S) 模 拟 月 径流 01 为 
224.31 ms ,下 热 面 变化 情景 (S;) 模 拟 月 径流 0; 为 
218.10 mx…s!', 协 同 变化 情景 (S;) 模 拟 月 径流 0; 为 
224.72 m…s"'。 人 气候 变化 条 件 下 (So 一 51) ,人 研究 区 月 
径流 量 增加 了 6.62 ms, 气候 变化 对 径流 影响 的 
贡献 率 (Be) 为 94.2%。 下 垫 面 变化 条 件 (S 一 S$,) 引 


表 6 气候 及 下 垫 面 变化 对 


起 的 径流 增加 量 为 0.41 m.s…, 下 垫 面 变化 对 径流 
影响 的 贡献 率 (BL) 为 5.8%。 在 气候 与 下 垫 面 变化 
的 协同 变化 条 件 下 (Ss 一 5;) ,研究 区 月 径流 量 增加 
7.03 ms*'。1968 一 2017 年 ,气候 变化 、 下 热 面 变 化 
及 二 者 协同 变化 均 引 起 研究 区 月 径流 量 的 增加 , 径 
流 对 气候 及 下 垫 面 变化 协同 作用 的 响应 具有 受 加 
效应 ,其 中 气候 变化 对 人 研究 区 径流 演变 的 影响 占据 
主导 地 位 。 


3 结论 与 讨论 


1968 一 2017 年 , 叶 尔 羌 河流 域 上 游 河源 区 呈 暧 
湿 化 的 气候 特征 ,加 剧 了 冰川 消融 的 速度 ,直接 导 
致 了 当地 径流 量 的 增加 。 本 文 对 叶 尔 羌 河 流域 上 
游 河源 区 气象 水 文 变化 特征 进行 了 分 析 ,基于 
SWAT RAI X} 1968—2017 年 的 月 径流 进行 模拟 ,区 
分 并 量化 了 气候 与 下 垫 面 变化 对 径流 变化 的 影响 ， 
主要 结论 如 下 : 

(1) 研究 区 近 50 a 的 降水 量 、 气 温 及 径流 量 增 
加 趋势 显著 ,年 均 增 长 率 分 别 为 1.0% 、1.2% 和 
0.5%。 人 研究 区 内 冰川 覆盖 率 达 12.6% , 除 降雨 和 积 
雪 外 ,在 模型 的 模拟 过 程 中 应 充分 考虑 冰川 水 文 过 
程 , 故 本 文 在 SWAT 模 型 中 乱入 了 基于 增强 温度 指 
数 法 的 冰川 模块 。 增 强 温度 指数 法 在 数据 缺乏 的 
流域 尺度 计算 上 具有 较 强 的 鲁 棒 性 , 且 该 算法 整体 
性 能 好 , 相 比 度 日 因子 法 ,该 法 能 更 好 地 解释 由 局 
部 地 形 引 起 地 空间 蜡 质 性 中 。 舱 入 冰川 模块 的 
SWAT 月 径流 模型 在 校准 期 .验证 期 与 全 时 段 的 模 
拟 结果 评价 均 为 “ 优 ”, 表 明 模 型 在 该 流域 具有 较 好 
的 适用 性 ,但 依然 与 实际 径流 变化 过 程 存在 一 定 差 
中 ,导致 这 一 差距 的 原因 包括 以 下 几 个 方面 :个 别 
年 份 存在 冰 湖 省 决 洪水 现象 2 ;模型 数据 在 输入 过 
程 中 存在 不 确定 性 以 及 数据 获取 具有 一 定局 限 性 。 

(2) 经 模型 模拟 得 到 的 年 际 间 的 年 径流 量 及 冰 
川 径 流量 结果 均 呈 增加 趋势 ,其 中 多 年 冰川 径流 量 
占 径流 总 量 的 51.19 ,该 结果 与 其 他 学 者 对 其 他 类 


径流 影响 定量 分 析 结 果 


Tab.6 Quantitative analysis results of the influence of climate and underlying surface change on runoff 


水 文 站 Own s" On s” Qm s" Q/M? s” 


AQcdm s 


AQUm':s" AQca/m s! Be Bd% 


卡 群 217.67 224.31 218.10 224.72 


6.62 0.41 7.03 94.2 5.8 


注 :0,、Q01、Q0..0; 分 别 表示 基础 情景 气候 变化 情景 下 垫 面 变化 情景 .协同 变化 情景 下 的 模拟 月 径流 ;AQc、AQL、AQc 分 别 表示 径流 在 气候 、 下 


垫 面 \, 气 候 和 下 执 面 共同 作用 下 的 变化 量 ;B. 和 B. 表 示 气 候 变 化 和 下 垫 面 变 化 对 径流 影响 的 贡献 率 。 
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似 气 候 条 件 下 的 冰川 融雪 补给 型 河流 的 研究 结 
相似 中 。 这 说 明 暧 湿 化 的 气候 特征 对 冰川 这 类 周 
态 水 资源 系统 的 影响 十 分 剧烈 , 且 随 着 暖 湿 化 特征 
的 不 断 显现 ,未 来 冰川 径流 在 总 径流 中 的 占 比 将 会 
逐渐 增 大 。 从 年 内 模拟 结果 看 ,月 冰川 径流 量 占 月 
径流 量 的 比重 在 11.9%~74.0% 之 间 , 其 中 2 月 最 小 ， 
8 月 最 大 ;研究 区 的 冰川 产 流 主要 发 生 在 6 一 9 月 , 占 
全 年 冰川 径流 量 的 90.0% 以 上 。 值 得 注意 的 是 ,50 a 
间 研 究 区 径流 量 的 峰值 期 提前 了 近 1 个 月 左右 (从 8 
月 提前 到 7 月 ) ,表明 当地 的 冰川 消融 速度 正在 不 断 
Apos 

(3) 气候 变化 和 下 垫 面 变化 分 别 导 致 研究 区 月 
径流 增加 6.62 ms 和 0.41 m3…s ,二 者 对 径流 变化 
的 贡献 率 分 别 为 94.17% 和 5.83%, 从 气候 和 下 热 面 
变化 对 当地 径流 影响 的 分 析 结 果 来 看 ,在 气温 升 高 
和 降水 量 增加 等 气候 条 件 变 化 的 影响 下 ,人 研究 区 模 
拟 月 径流 量 有 较 大 增加 ;在 下 垫 面 变化 的 情景 下 ， 
模拟 月 径流 的 增加 量 相对 较 少 ,可 见 气 候 变化 对 研 
究 区 径流 变化 的 影响 占 主导 地 位 。 

需要 说 明 的 是 ,本 研究 还 存在 一 定 的 不 足 :(1) 
本 研究 区 的 气象 数据 集 获取 具有 一 定 难 度 , 今 后 考 
虑 收集 基于 遥感 资料 同化 的 再 分 析 数 据 集 并 应 用 
于 后 续 模 拟 研究 。(2) 未 考虑 到 水 利 工 程 “ 等 因素 
对 人 研究 区 径流 的 影响 ,在 日 后 的 模拟 过 程 中 应 对 
这 些 人 为 因素 的 影响 程度 进行 更 深入 的 人 研究”。 
(3) 未 结合 未 来 变化 情景 下 不 同 的 气候 模式 ,预测 
研究 区 径流 的 演变 过 程 ,在 日 后 的 研究 中 应 进一步 
对 气候 变化 条 件 下 径流 的 变化 规律 进行 分 析 ,为 流 
域 未 来 的 水 资源 管理 提供 决策 支撑 。 
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Abstract: Climate and underlying surface changes are the main factors affecting the changes in the hydrological 
cycle of a river basin. Research on the response of the runoff evolution to climate and underlying surface changes 
at different temporal and spatial scales is the basis of water resource planning and management. This study took 
the upper reaches of the Yarkand River Basin in Xinjiang, China as the study area, constructed a SWAT 
distributed hydrological model embedded with glacier modules, and simulated the monthly runoff from 1968 to 
2017 at the Kaqun Hydrological Station to evaluate the model' s applicability in the study area. The simulation 
results were used to analyze the interannual variability and the intra-year distribution of the glacier runoff in the 
study area and quantitatively calculate the contribution rate of climate and underlying surface changes to runoff 
changes. The results show that the model has good applicability during the runoff simulation in the study area. 
The R° during the calibration period from 1968 to 1992 and the verification period from 1993 to 2017 was 0.77 
and 0.86, respectively; the NSE was 0.76 and 0.85, respectively; the RSR was 0.49 and 0.38, respectively; and 
the PBIAS was -7.4% and 0.6%, respectively. The simulation showed that the glacier runoff in the study area 
generally increased in the past 50 years, with the proportion of the total runoff being approximately 51.196. The 
glacier runoff mainly occurred from June to September of the year, accounting for 90.095 of the annual glacier 
runoff. The changes in climate and land use increased the monthly runoff by 6.62 m°- s ' and 0.41 m°- s”, respectively. 
The contribution rate of climate change is 94.296, depicting that the impact of climate change on the runoff in the 
study area is dominant. The research results can provide a scientific basis for understanding the causes of the 
historical runoff changes in the study area and predicting the future runoff evolution trends. 
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